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Abstract. We consider perturbative modifications of the Friedmann equations in terms of energy 
density corresponding to modified theories of gravity proposed as an alternative route to comply 
, with the observed accelerated expansion of the universe. Assuming that the present matter content 

■ of the universe is a pressureless fluid, the possible singularities that may arise as the final state of the 
universe are surveyed. It is shown that, at most, two coefficients of the perturbative expansion of the 
Friedman equations are relevant for the analysis. Some examples of application of the perturbative 
scheme are included. 
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*> INTRODUCTION 

in 

Recent astronomical data from Type la supernovae 0]] as well as from the CMB spec- 
■^j- ■ trum [0] confirm that our universe is undergoing an accelerated expansion period. In 

^ . order to comply with this feature, one may resort to postulating a dark energy content 

for the universe yQ, with undesired properties, such as violation of some energy con- 
ditions, or going beyond general relativity in the quest for another theory of gravity 

With either approaches for dealing with the observed accelerated expansion, cosmol- 
ogy is much richer than it was thought in the previous century. According to classi- 
cal cosmologies, the universe started at an initial singularity, the Big Bang, and it was 
doomed to expand forever, since the matter content was not dense enough to stop expan- 
sion and collapse into a final singularity. 

However, accelerated expansion of the universe has lead to consider other plans for 
the end of the universe in the form of Big Rip singularities [8], directional singularities 
[@] or milder sudden singularities lHoll . It is therefore interesting to tell under which 
circumstances a modified theory of gravity may lead to such fate for the universe. 



MODIFIED GRAVITY 

Instead of dealing with a full theory of modified gravity, we focus on the consequences 
at the level of cosmological equations, namely Friedmann equation, just requiring that it 



admits a generalised power expansion on the density p around a value p*, 



(?) = H 2 = h (p-p^ + h l (p-p^ + ---, &<^i<.... (1) 

The standard density term arises as the lineal term with an exponent equal to one and 
the cosmological constant appears with a null exponent in this expansion. Further terms 
are interpreted as modifications of the theory. 

On the other hand, the energy conservation law implies 

p+3H(p+p)=0, (2) 

but assuming that the accelerated expansion of the universe is solely due to the modi- 
fication of the theory of gravity, we choose a pressureless dust as matter content of the 
universe, so that the scale factor of the universe and the density are related by pa 3 = K. 

Thereby we may get rid of the scale factor and write down a modified Friedmann 
equation in terms of density: 

£ = -3^(p-p*)^ 2 -|^(p-p*) M/2 + - • O) 
p 2 y/ho 

Solving this equation provides a perturbative expansion of the density in coordinate 
time, which we want to compare with a similar expansion for the scale factor: 

a(t) = c \t-t \ 7l0 + ci\t~t l \ 71] +■■■, TJo < TJi < 

in terms of coordinate time. This is useful, since in iflTll we have related the exponents 
r\i with the strength of singularities according to the standard definitions by Tipler [Q3] 
and Krolak II 1311 . as we see in Table 1. The strength of the singularities just points out if 
tidal forces are strong enough to disrupt finite objects on crossing them. 

TABLE 1. Singularities in cosmological models 
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The last column refers to the classification of future singularities in [Q3]): 

• Type I: "Big Rip": divergent a. 

• Type II: "Sudden": finite a, H, divergent H. 

• Type III: "Big Freeze": finite a, divergent//. 

• Type IV: "Big Brake": finite a, H, H, but divergent higher derivatives. 



We have chosen the decreasing density branch of the Friedmann equation, since we 
wish to mimic the actual expansion phase of the universe. There are two possibilities: 
either there is a critical value p* for which the density stops decreasing or it goes on 
decreasing without a lower bound. The results are consigned in tables 2 and 3. More 
details may be found in lfl5ll . 

Singularities in models with null critical density are determined by the sign of the 
first exponent If it is positive, expansion goes on forever and no singularity arises, 
but if the exponent is negative, a Big Rip singularity comes up due to the accelerated 
expansion. 

TABLE 2. Singularities in models with p* = 
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On the contrary, the structure of models with a finite critical density p* is much 
richer. No Big Rip singularities appear and they are all weak. Therefore, they cannot 
be considered the final stage of the universe iflitl . 

Models with negative £o exponent have divergent H, which is a Big Freeze singularity, 
which is weak under Tipler's definition, but strong with Krolak's. 

Models with null £o exhibit a cosmological constant term and produce just weak 
singularities. The second exponent £\ can be used to tell the derivative of H which is 
singular, but none of these can be taken as the final fate of the universe, due to the 
weakness of the singularity. 

Finally, models with no cosmological constant, £o > show the same types of weak 
singularities. 



Examples of proposed models belonging to these families may be found in II 1511 . 
TABLE 3. Singularities in models with p* 7^ 
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